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We investigate the fragmented many-body ground states of a spin-2 Bose gas in zero magnetic field. We point 
out that the exact ground state is not simply an average over rotationally-invariant mean-field states, in contrast 
to the spin-1 case with even number of particles N. We construct the exact ground states and compare them 
with the angular-averaged polar and cyclic states. The angular-averaged polar states fail to retrieve the exact 
eigenstate at W > 6 while angular-averaged cyclic states sustain only for N with a multiple of 3. We calculate 
the density matrices and two-particle density matrices to show how deviant the angular-averaged state is from 
the exact one. 
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Introduction.- Since the advancement of optica^ trapped 
Bose-Einstein condensate (BEC) lU], spinor BEC Q has pro¬ 
vided a paradigm to study magnetism, spin textures, topolog¬ 
ical excitations, and quantum dynamics of associated many- 
body ground states The quantum phases of spin-/ BEC 
can be ferromagnetic or polar for / = 1, 2 or cyclic {f= 
2)|l7H depending on two-body s-wave scattering lengths ap 
of even total spin F' up to 2/. Uniaxial and biaxial spin ne¬ 
matic phases can also be identified in spin-2 case when the 
degeneracy in polar phase is lifted by thermal or quantum 
fluctuations lUOlflin . Higher spin Bose gas can involve even 
more complicated phases |[l2l - ll^ . Rich spin mixing dynam¬ 
ics has been used to observe the ferromagnetic |[l^ and anti¬ 
ferromagnetic (polar) properties in respectively for spin-1 
^'"Rb and ^^Na Bose gases. In spin-2 cases, polar phase is 
the likely phase for ^'’Rb OlM [l^ . Recent experimental de¬ 
velopments in^inor BEC involve spin textures m . spin 
dynamics 11^ m under quadratic Zeeman shift. There are 
also studies of quantum phase transitions by Earaday rotation 
spectroscopy B or adiabatic microwave fields B . and spin 
coherence measurements by Ramsey interferometry 11^ 1^ . 

A condensate of bosons forms when one of its single¬ 
particle wavefunction is macroscopically occupied Iz^ . The 
fragmentation of BEC becomes feasible when multiple 
macroscopic single-particle densities are degenerate in spinor 
Bose gases 1^ though it is fragile in the presence of weak 
external magnetic fields or symmetry-breaking perturbations 
1^ . Lately many interests in fragmented BEC include dy¬ 
namical formation of two-dimensional fragmented BEC 11^ . 
quadratic Zeeman effect on spin fragmentation and 

fragmented many-body ground states with anisotropic long- 
range interactions 1^ or trapping potentials It has been 
proposed that signatures of fragmentation can be probed by 
measuring density-density correlations 1^ while fragmen¬ 
tation resulting from Goldstone magnon instability 11^ and 
spin-orbit coupling |j^ are also investigated. 

The frag mented structure of the ground state in spin-1 Bose 
gases originates from the rotational invariance 

in spin degrees of freedom. The symmetry-breaking mean- 
field (ME) treatment fails in describing the exact ground state 
for presumption of single spin coherent condensate. Eor scat¬ 
tering lengths obeying 02 > oq, the ME state is polar, but 


the exact ground state is fragmented and, for even number 
of particles N, can be viewed as a collection of two-particle 
spin-singlets. This exact ground state has equal populations 
in the magnetic sublevels with large number fluctuations of 
order N pM . which is very different from ME states. It is 
claimed that this exact ground state can be understood 
as the angular-average of the ME polar states as an analog 
to the relation between Lock and coherent states in a double¬ 
well system 1^ . This remains the view adopted by the most 
recent review articles Hi. Is this perspective of angular- 
averaged states universal and applicable in constructing the 
exact ground states for larger spins? In this paper, we inves¬ 
tigate many-body ground states of a spin-2 Bose gas, and 
demonstrate how angular-averaged states are unable to con¬ 
struct them. That the angular averaged ME states is the exact 
ground state is just a coincidence in spin-1 system. We ad¬ 
dress the inapplicability of the angular-averaging process, and 
also show how the angular-averaged ME state deviates from 
the exact eigenstates by studying the two-particle density ma¬ 
trices. 

Spin-2 Bose gas.- Eor a spin-/ Bose gas at low tempera¬ 
ture, the two-body particle interaction involves only scatter- 
ing channels of even total hyperflne spin F' states up to 2/ 
ra . We shall consider the single-mode approximation (SMA) 
where the spatial part of the wavefunction is the same for all 
spin sublevels such that the field operator ipm (r) = \/ p{r)dm 
with the density p{r) and spinor operator dm- Since the spa¬ 
tial part is frozen, the effective Hamiltonian (in zero magnetic 
field, to which we shall limit ourselves) involves only the in¬ 
teraction V which reads H 


V = 
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FIG. 1. (Color online) Mean-field-plus (MF’'") phase diagram and many-body ground states of spin-2 Bose gas in parameter space ( 7 , /3). (a) 
Mean-field-plus phase diagram; ferromagnetic (F), polar (PO and P2), and cyclic (C) phases, (b-i): phase diagram for exact many-body states 
with finite number of particles from N = 2 — 9. The form of the many-body ground states are shown. These states are constructed by spin- 
singlet pairs (©2) and trios (©3), and in certain regions, creation operator oj (2I22) for a single particle (pair of particles) with / = 2 , m = 2 
are also necessary (see text for definition of these operators). States that involve a\ or explicitly are degenerate with their partners obtained 
by rotational symmetry (not shown). Dots in the formulas, when the complete expressions are not given, indicate that linear superposition with 
other terms are required. Wavefunctions shown here are not normalized. Dashed lines indicate the (schematic) phase boundaries. 


where the coefficients are ct = ( 432 + 354 )/?,/3 = ( 34 — 32 )/?, 
and 7 = (30 — 34)/5 — 2(32 — 34 )/?. Here the interaction pa¬ 
rameters 3 f = 4,TTh^aF/M with the mass of the atom M and 
s-wave scattering length ap, and (00|2n; 2n') is the Clebsh- 
Gordan coefficient for the overlap between the states with two 
spin -2 bosons of = n, n' and the spin singlet | 00 ). 

In MF theory, bosons condense. Particles macroscopically 
occupy a single quantum state which can be described by a 
spin-2 wavefunction (v?- 2 , ■■■, ¥^ 2 ). In our case, there are three 
phases characterized by two order parameters of magnetiza¬ 
tion (/) = Y^m =-2 and spin-singlet pair amplitude 

( 02 ) §1 , where 

2 

02 = ^ V5{00\2m;2 — m)ama-m 

m——2 

= 2 ( 120,-2 — 2aia_i + Og (2) 

is an operator which annihilates a singlet pair. That is, we 
have (© 2 ) = Thsre are three phases. 

The ferromagnetic (F) phase has a finite (/) and zero (© 2 ), 
while the polar phase (P) has (© 2 ) = 1 without (/). When 
/3, 7 > 0, the cyclic (C) phase has the lowest mean-field en¬ 
ergy for both zero (© 2 ) and (/), breaking the time-reversal 
symmetry. The phase boundary between F and polar phases 
is delineated by the line 4/3 = 7. The phase diagram is as 
shown in Fig. [Ha) flUtl- Representative wavefunctions are 
( 1 , 0,0, 0,0) for F; ( 1 , 0,0, 1) which is equivalent by ro¬ 

tation to (1,0, v/2, 0,-1) for C 13 • Within mean-field, the 


polar phase P can have wavefunctions PO = (0, 0,1, 0,0) orP2 
= (1, 0 , 0 , 0 ,1 ), or any real linear combinations thereof 
(apart from rotations). This degeneracy however is lifted by 
fluctuations id [3. The resulting phase diagram, which we 
shall call the mean-field-plus (MF+) phase diagram, is shown 
in FigllJa). We note that P0(2) can be also represented by the 
polynomial forms via spherical harmonics as ( 2 ^^ — — 3 ^) 

and (x^ — 3 ^) respectively 13 dl- For ease of referral later, 
we shall call the regions in (/3, 7 ) parameter space occupied 
by the F, C, P phases as F, C, P regions. 

Many-body ground states.- Let us now discuss the many- 
body ground states of a spin-2 Bose gas in zero magnetic field 
with SMA. ||8l|3,|3- For a given N, the many-body ground 
states are characterized by two quantum numbers F and r. 
F is the total spin and the integer quantum number r can be 
interpreted as the number of particles other than spin-singlet 
pairs, therefore r is given by Srigo + A 11113 where ngg is the 
number of spin-singlet trios and the integer A indicates parti¬ 
cles other than spin-singlet pairs and trios. These states are 
also eigenstates of the operator A = 0|02 with eigenvalues 
A = N{N -F 3) — t(t + 3) = {N — t){N + 3 + t), from 
which T can be evaluated. The exact ground state energy is 
proportional to PF{F + 1) — 7 t(t -F 3) aside from a term 
depending only on N. The phase diagram derived by mini¬ 
mizing the ground state energy is as sketched in Fig. [T]for N 
= 2 to 9. The line that separates the two phases in the /3 < 0, 
7 < 0 region is given by {4N -F 2)/3/(W + 3) = 7 13 which 
approaches to MF phase boundary in thermodynamic limit. 

The wavefunctions listed in Fig [T] are constructed accord- 
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ing to the eigenvalues F and r which minimize the energy. 
In the ferromagnetic region F, they are (d 2 )^| 0 ), hence iden¬ 
tical with the mean-field states. (Here |0) denotes the vac¬ 
uum). For the polar P region, wavefunctions differ according 
to whether N is even or odd. For even N, the wavefunctions 
are ( 02 )'^'^^|O) corresponding to r = 0 , F = 0 , and maximum 
possible A’s which are N {N + 3) (see also iisll ). For odd N, 
the states that appear near the —7 axis have wavefunctions 
(up to rotations) again correspond to states 
with maximum possible A’s which are now (TV — 1)(TV + 4 ) 
with r = 1 and F = 2. Near the +P axis but still 7 < 0, the 
states (for N > 3) have the form 03 ( 02 )^'^“^^'^^|O) where 

03 = -VQalalal^ + + |(a-i)^4 

- (3) 

is an operator which creates a spin -0 trio (not normalized). 
These states have r = 3 (see the end of Sec IV in ll43l] ) and A 
= {N- 3){N + 6 ) with F = 0. 

In the region C, if the particle numbers are multiple of 3, 
the ground state wavefunctions are singlets constructed by ©J 
and ©2 with the eigenvalue A = 0. The wavefunctions are 
indicated in Fig [l1 an d their derivation can be found in Supple¬ 
mental Material For other N"s, the C region is divided 
into two parts. The states near the +/3 axis and 7 > 0 are 
again spin singlets. If = 2(mod3), hence N = 3R + 2, 
they are given by ©JIIPsh) (with r = 3R and F = 0) where 
is the corresponding ground state in region C for 3R 
particles. For N = l(mod3) and hence N = 3R + 4 for inte¬ 
ger R (when TV > 4 ) then the wavefunctions are (02)^|'I'3fl) 
and again r = 3R. For the TV’s shown in Fig[T] these states 
happen to be the same as the states near the +/3 axis on the 
7 < 0 side so that they are the same phase, but this needs not 
hold for larger particle numbers TV > 11. The states near the 
+7 axis with ^ > 0 have instead finite magnetization and are 
not rotationally invariant. They are not directly relevant in the 
rest of the paper and we shall not discuss them in detail. 

Angular-averaged mean-field states.- Now we turn to the 
angular-averagedMF states and compare with the correspond¬ 
ing exact eigenstates. First we recall the corresponding results 
for spin-1. The MF ferromagnetic state corresponds to the ex¬ 
act solution. Their angular average actually vanishes. The po¬ 
lar mean-field state has a finite average only for TV even, and 
give the correct exact many-body state . We then demon¬ 
strate how for spin -2 the averaging process enables the frag¬ 
mentation in both the polar and cyclic phases but in general it 
fails to correctly construct the corresponding exact eigenstates 
at a given point (/3, 7 ) in parameter space. 

The situation for the ferromagnetic state is exactly the 
same as the spin-1 case. We now consider the angular- 
averaged polar state of P0(2). Starting from the reference state 
(0,0,1,0, 0), the state obtained by rotations via the Euler an¬ 
gels a, /3, 7 , which we denote collectively as H, is given by 
where the matrix is the spin- 

2 irreducible representation of the rotation operator 114^ (see 


also Supplemental Materials ll43l] ). The general (unnormal¬ 
ized) rotationally invariant state is constructed via 

= im. /ft*”’'"”" 

where we have defined ^ s Jf /» ^ /“ g. For 

our polar state PO, we thus use a’l'(H) —>• We 

call the resulting state |4'po)ai;- It is straight-forward to eval¬ 
uate the angular integrals. We find that |t['po)ai; retrieves the 
exact ground states for even TV = 2,4. For odd N = 3,5, 7, 
we recover the spin singlet state located near the -\-f3 axis as 
indicated in Fig. [T](see more details in ll43l] ). (obviously the 
angular average cannot produce the states with finite magne¬ 
tization near the —7 axis). However for even TV > 6 and odd 
TV > 9, the angular-averaged states fail to construct the exact 
ground states. For example, for TV = 6 it gives rather 
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which is in fact not even an eigenstate of the Hamiltonian in 
Eq. O. 

Similarly for P2, we can construct angular-averaged states 
as in Eq. © except now we use dl( 0 ) —^ 
with + ^m!- 2 )(^)- In this case the an¬ 

gular average vanishes if TV is odd. Eor TV even again it pro¬ 
duces the correct ground states for TV = 2,4 but fails again at 
6 . 

Actually why the angular averaged mean-field states can or 
cannot produce the manybody state is now clear. Eor TV up to 
5, the exact many-body singlet states are unique. Since angu¬ 
lar averaged mean-field states must either be zero or they must 
be a rotationally invariant, they must either vanish or produce 
the singlet states. This is actually independent of whether the 
starting mean-field state is the corresponding ground state for 
the given parameters in the Hamiltonian. Eor TV = 6 , 8 , 9.., 
the many-body singlet states are no longer unique. The angu¬ 
lar average, if it is not zero, just produces some linear combi¬ 
nations of these singlets. The resulting states have nothing to 
do with the ground state solutions of the Hamiltonian. That 
the angular average of the polar state for spin -1 produces cor¬ 
rectly the exact many-body state for even TV is purely because 
that, for spin- 1 , this singlet is unique. 

Let us also consider the angular-averaged states for a linear 
combination of both P0(2), and use = cos0D^\(Cl)-\- 

sin 0 (F^^ 2 (i^) F /\/2mEc{. @. Eor TV = 6 , and 

the angular-averaged polar state becomes (see El) 
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which in general again fails to become an eigenstate of Eq. 
O. The 9 dependence of the coefficients obtained above can 
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be understood by the symmetries of the general polar state un¬ 
der 0 —9 and 6 —>■ n/S — 6 ll43l] . We note that 

never produces the many-body state — i^( 02 )^]|O) 

in the C region. It happens that when cos 60 = —1/3, the 
angular-averaged polar state becomes the exact ground state 
(© 2 )^ 10 ) in the P region. 

The above special value of 9 can be understood as follows. 
It can be shown that the weighted average 3/2 d9 sin 30 
over 0, together with the average over Euler angles above, is 
equivalent to an average over the 4-sphere in the quantum ro¬ 
tor picture of If we apply this average to ip^{Cl), we 
obtain the exact many-body state |0) for even N (the 

average vanishes for odd N). This is because the above men¬ 
tioned averages guarantee that we obtain a state that is invari¬ 
ant under SO(5) rotations, and (©|)^/^|0) is the only such 
state (corresponding to r = 0 of 1421] '). The 0-averaged of 
cos(60) is —1/3. 

The comparison between the angular-averaged polar state 
and the exact eigenstates can also be viewed in a different 
manner. Let us consider the operator ©^© 2 . We note that, for 
the polar state, = 1 for any Ct and 

hence, for the normalized state |^p( 0 ))aD of |4'p(0))a,;, we 
have the expectation value 

a,(’l-p(0)|0|02|^p(0))a. = N{N - l)XN-2{e)/XN{9), 

( 6 ) 


where Xn{9) = a«(^'p(0)|'I'p(0))a« and we have defined 
Xq = 1. While this formula is general, let us focus on PO. 
The evaluation of at at 0 = 0 are particularly straightfor¬ 
ward. We have 


Xw(0)= a.(^p(0)|^'p(0))a„ 



where 11 = ^172 represents the rotation 112 followed by the 

/o'! '' 

inverse of Hi. Here we have used the relation Dq g(H) = 

S. We obtain X^ = 0, X 2 = 1/5, 
X 3 = 2/35, X 4 = 3/35, X 5 = 4/77, Xg = 53/(7 • 11 • 
13), X 7 = 6/(11 • 13), Xg = 5 • 19/(11 • 13 • 17), Xg = 
2^ • 197/(11 • 13 • 17 • 19). On the other hand, as already 
mentioned, the exact eigenstates are also eigenvectors of the 
operator A = ©^© 2 . For even N, the states ( 02 )^'^^|O) have 
eigenvalues A = N{N -f 3). For odd N, the exact eigenstates 
( 02 )^^“^^^^ 03 |O) have eigenvalues (X —3)(X 4 - 6 ). We can 
check directly from Eq. (| 6 ]) that the expectation values for 
©202 equal these exact values for N = 2,3, 4, 5, 7 but not 6 , 
8 or 9. 

Cyclic states.- We now study the angular-averaged cyclic 
states. It is simplest to use the reference state :^(1, 0, 0, 

^ 2 , 0 ) and hence (p^{n) = i" 

Eq. (0]). The angular-averaged C states are finite only when 
X is a multiple of 3, which can be easily seen by consider¬ 
ing the integral over the angle 7. It turns out that, in these 
cases, the angular-averaged C states do produce the correct 
many-body states. This is due to the fact that (H) obeys 




FIG. 2. (Color online) Two-particle density matrices for the angular- 
averaged polar state for X = 6. Three density matrices for J = 0, 2, 
4 are shown in (a), (b), and (c) respectively. 

~ hence the angular aver¬ 

aged state is annihilated by © 2 , so that the resulting state must 
satisfy A = 0 and hence correctly produce the corresponding 
many-body state. We have also verified this conclusion by 
direct angular averages (see S). 

Reduced density matrix.- As a further investigation, com¬ 
pare, for X = 6, the two-particle density matrices for the 
angular-averaged polar states with those for the exact many- 

body states IT's) = — 7 ^=(© 2 )^| 0 ) where ^(3) is a normal- 
s/ais) 

ization constant. (The one-particle density matrices are obvi¬ 
ously identical since both states are rotational invariant). It is 
simplest to present the results using the operators 

Ajm = ^ ^ {^J^d[\22TTi\Tn‘2^(Xjyi.^(Xjyi2 • (8) 

mi ,7712 

(T'gjAj^Aj/M'l^e) is finite only when J = J'and M = 
M', and is further M independent, as expected by rotational 
invariance. These values are discussed in 1^ . We have 
(T' 6 |i^;o^o|^ 6 ) = 54/5 = 10.8, and = 

(4'6|ilMA4M|«'6) = 48/35 1.37. 

The numerical results for the angular-averaged ME state is 
shown in Fig|2l The values oscillates with 0 with period 7 r /3 
due to the cos 60 factor in Eq. (|5]). For general 0, the differ¬ 
ence between the angular-averaged ME and the spin-singlet 
pair states is less than 10%. For example (A^gAgg) = 10.8 
for the exact many-body state while (AjgAgg) = 9.7 in Fig. 
|2a) at 0 = 0. The values are identical at cos(60) = —1/3. 

While the density matrices at finite X in general differ, 
it can be shown 11 ^ that they have the same leading terms 
in the large X limit, so that the energy per particle remains 
the same up to corrections of order 1/X, as in the case for 
spin-1 ll2^ . In the large N limit, the fragmented state has 
macroscopic number fluctuations while they decay rapidly as 
miniscule magnetization sets in, therefore it is fragile against 
symmetry-breaking perturbations. However we expect that the 
fragmentation of many-body ground state can be observable 
in the few-particle system where its signature of two-particle 
correlations is more noticeable in contrast to the mean-field 
results. 
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In conclusion, the many-body ground states of spin-2 Bose 
gas in zero magnetic field are in general fragmented which is 
however not describable via angular-averaged MF states. For 
polar states the angular-averaged calculation fails to describe 
the exact eigenstates when even or odd iV > 6 or 9. For cyclic 
states, the angular-averaged treatment only sustains the exact 
ground states for particle number of a multiple of 3, which 
preserves the constraint of (©2 02) = 0. That the angular- 


averaged MF states for even N in spin-1 Bose gas are equiv¬ 
alent to the exact ground states is simply a coincidence. For 
even higher spinor BEC (/ > 3), we expect angular-averaged 
states fail to retrieve the exact eigenstates at even smaller num¬ 
ber of particles. 
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SUPPLEMENTAL MATERIAL OE ERAGMENTED 
MANY-BODY STATES OE SPIN-2 BOSE GAS 
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In this section, we reproduce the spin-2 irreducible 
representation of the rotation operator D^} = 

a+m 7 )^^) jidl] for Euler angles a, /3, 7 . In ma¬ 
trix form, is 




m' ,m \ 

/ 


cos^d) 


sin /3 


(1 -I- cos /3) 


|sin^/3 ^d^(cos/3 —1) 


sm 


4 /3 


-I-cos/3) |(2cos/3 — l)(cos/3-f 1) —-i /1 sin/? cos/3 |(2 cos/3-I-1)(1 — cos/3) ^d^(cos/3—1) 


I sin^ /3 


I sin /3 cos /3 


^(3cos^ /3 — 1) 


— W I sin/3 cos/3 


I sin^ /3 


—cos/3) i(2cos/3-f 1)(1 — cos/3) ^1 sin/3 cos/3 i(2 cos/3 — l)(cos/3-I-1) — ^^^(1-I-cos/3) 


V 


sm 


4 /3 


sin /3 


(1 — cos^) 


sin2/3 ™^(cos/3 + l) 


COS 


4//3 


, (9) 


(f) y 


I- 

which is expressed in terms of spin bases {ip 2 , ipi, (fio, the wavefunction becomes 


I. ANGULAR AVERAGE MEAN-EIELD POLAR STATES 
FOR FINITE NUMBER OF PARTICLES 


From Eq. (4) in the paper and with a^(ii) = 
where = cos 6 »£)^U^) -f 


sme{Dl^] 2 {n) + dI^]_ 2 {Ci))/V 2 , we first 
a and 7 , which gives 




N 


^ / N'\ cos" 0 sin ^ 9 /N — n 

= ( „ ) 2(^-")/2 I ^ 

n—0 ^ z 






( 11 ) 


where (iV — n)/2 is integer. We then evaluate the /3 average 
average over either analytically or with the help of Mathematic a. In the 
below, we report the results for this angular-averaged polar 
states for finite number of particles At = 2 to 10 . 


\'^'p{0,P))o 


N 


= E 

n—Q 


N' 




N-r, 


N-n 


X \y jp)a\ +y ,(/3)ad 

I / ^ 7712,2''^ ' m 2 / ^ rrig,—2^^ ‘ m^ 


7 ( 2 ) 

^ ^/(m 7 , 0 ^/(m ),0 


( 10 ) 


where /(m) = J2m=mi 2 expand the above fur¬ 

ther and use one of the delta function constraint (3/(m),o> 


A. N=2 

From Eq. 0 , we have the angular-averaged polar state 
1 


\^p{e))^^ = y= I \%i9,/3)}smpdp, 

25.20-2 ~ SfijaLi + (oq)^ | 0 )) 


2V^.Jo 

1 


5 V 2 ( 


( 12 ) 


Note that it has no 9 dependence. The angular averaged ME 
state reproduces the exact many-body state. 
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B. N=3 


E. N=6 


From Eq. (HB, we have the angular averaged polar state 

l'I^p(^))a. = ^ £ sinPd(3, 

= ^ (cos^ 0 — 3 cos 9 sin^ 6 ) © 3 10) 

OO 

= A[cos(30)]0j|O), (13) 

where ©3 is a three-particle singlet operator, 


For this even number of particles, we expect a combination 
of two- and three-particle singlet states to appear. From Eq. 
(HB, we have 


l'I'pW)a„ = 


1 


2 V&., 

1 


lOOlv^ . 

-f (12-f 36 cos 60) ( 0 ;^) 


sin/3d/3, 

(47- 2 cos 60) (©1)^ 
2 ' 


| 0 ). 


(17) 


©3 = ■^*^“ 0 )^ - -^alaja-i + + ^4(a-i)^ 


-=aoaAa 




2U.Qii_2. 


(14) 


Note that this angular averaged state has 0 dependence with 
a period of 27r/3 but always reproduces the exact many-body 
state for the region 7 < 0 and near the +(3 axis. 


We may express this wavefunction in terms of the normal¬ 
ized many-body state of two- and three-particle singlet states 
for N particles. 


(2) \ 
N=6 / 

- . ^_ f02)^O)> 

V24.33-5-7 V V 

00 

(3) \ 
N=6 / 


(19) 


For 0 = 0, the normalized angular-averaged state is 


C. N=4 

For this even number of particles, we again have 

where so far we still have a 0 independent angular average. 
This state which is a N/2 spin-singlet-pairs state, i.e, the exact 
eigenstate. 


^p(O) 


av 

1 


Vll •13-53 L 


5 • 3®/^ 


4- 


( 2 ) 

tV=6 


-^2^ • 7 




(3) 

Af=6 


( 20 ) 


where we note the finite overlap = 2/{\/3-7). 

The angular average has 0 dependence in general, and it can 
be expressed in terms of two-particle-singlets state only when 
cos 60 = —1/3. In general the angular-averaged polar state 
fails to construct the exact ground states which should be iV/2 
spin-singlet-pairs state. 


F. N=8 


D. N=5 


From Eq. (fTB . only n=l, 3, 5 are possible for (iV — n)/2 
is an integer. The angular averaged polar state is 

= ^£\%i9,/3))sinPdp, 

= —(cos® 0 — 3 cos 0 sin^ 0 

77^ 

— 2 cos® 0 sin^ ^)© 2 ® 3 | 0 ) 

= [cos(30)] ©^©I^IO). (16) 

Note that this angular averaged state again has 0 dependence 
with a period of 27r/3 and reproduces the exact many-body 
state for the region 7 < 0 and near the +13 axis. 


To investigate the 0 dependence of even number of parti¬ 
cles, we proceed to calculate the angular averaged polar state 
ofiV= 8 , 




1 


2\/^, 

1 


'^p{S,P)) sin f3d/3, 

(71 -8 cos 60) ( 0 ^)"^ 
2 ' 


243lV^ . 

16(3-f 9 cos 60)0^ (0;^) 


| 0 ). ( 21 ) 


Using the normalized singlet states. 


(2) \ 
N=8/ 

- . ^ _ (©4^0), 

V27.33-5-7-11 V V 

(22) 

(2.3) \ 
N=8/ 

- ,_^ 

V2-5-7-11-79 V 3/ 

(23) 
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we may express the normalized angular-averaged state (con¬ 
sider 0 = 0) as 


^p(O)) = 


3v^ L 


V3 




( 2 ) 

iV=8 


V79 




( 2 , 3 ) 

Ar=8 


(24) 


where again we use (V3 • 79). Note that 

when cos 66 = —1/3, the angular averaged MF state becomes 
the exact many-body state. 


G. N=10 


We may further investigate the angular averaged MF state 
for even N. From Eq. (dUl, we have 


I^PW)a. = 


2^10!. 


|^p(^')/3))sin/3d/3, 


(X 


(101-20 cos 60) 
-f 120(1-f 3 cos 60) 


|0). (25) 


Note that a special angle of cos60=—1/3 appears similar to 
the cases of N=6, 8. 

The 0 dependence obtained above can be understood as fol¬ 
lows. In the Cartesian representation, the general polar state 
is cos d{2z^ — — y^)j -b sin 0(a;^ — y^) 1^/2. 0 —>• —0 is 

equivalent to interchanging x and y, whereas 0 —^ tt 13 — 6 
has the effect of interchanging y and z as well as a sign 
change in the wavefunction. It follows that the angular av¬ 
eraged state must be invariant under 0^—6, while under 
0 —7r/3 — 0, it is multiplied by (—1)^. On the other hand, 
for N particles, the 0 dependence comes from terms of the 
form cos" 0 sin^“" 0 where n = 0, ...N with real coefficients 
(see Eq. (fTTI) ). Hence it must be of the form J2k=-N 
where c-k = cl.. Eor even TV, it follows that there is no 0 de¬ 
pendence for < 4, and the 0 dependence for 6 < N < 10 
can only be a linear combination of a constant and another 
term oc cos(60) (only c±6 and cq are allowed). For odd N 
with 3 < N < 7, the 0 dependence is via cos(30) (only c ±3 
allowed). 


II. ANGULAR AVERAGED MEAN-FIELD CYCLIC 
STATES FOR FINITE N 


When we angular averaged the mean-field cyclic state, we 
obtain, for N = 3, 

\'ifc)av = ^^:^03|O) (26) 

For TV = 6, we have 

" 7-11.13V^ 

(27) 


{qIY+ 18 [eiy 


| 0 ), 


In both cases, we produce the exact many-body states in the C 
region. 


III. WAVEFUNCTIONS FOR THE SINGLET MANY-BODY 
STATE WITH A = 0 IN THE C REGION 

We show how to obtain the singlet wavefunctions in the 
C regions of Fig 1 in the main text. To simplify notations, we 
shall often simply write (00|n, n') for the Clebsch-Gordan co¬ 
efficients (00|2n; 2n') when no confusion arises. We observe 
that 

m—— 2 

2 

= V5c {00\-m,m)ai^Alm, (28) 

m—— 2 
1 / 2 

where c = — ^ (|) . Such a relation is expected since both 

sides create a singlet state of three particles. It is easy to see 
that 

[dm, ^2 m] ~ 2{2M\m, M — 

= 2(-l)™(2TW-m|-m,TW)aE„. (29) 

It is useful to note that though [dm, ^ 2 m] d, we have 
2 

E[dm,iL]=0. (30) 

m—— 2 

This relation is expected since the left hand side is ro- 
tationally invariant but its right hand side can only in¬ 
volve one creation operator. Indeed, J2m=-2l^rn, = 

2 ^^^_ 2 (—l)^(20|bt, —y)o,Q but the sum is proportional to 
^^^_ 2 ( 00 |p, —p)(20|/i, — /r) = 0 due to the orthogonality 
between the states jOO) and |20). From Eq. ( l28l l we can eval¬ 
uate 


[d_m,0^] = -(-!)' 


. (3x7)V^ ^-t 


2,7715 


(31) 


and hence 


[02,0i]=-^^^^ E {^L,dm}, 

771 —— 2 
2 

=-(3 X 7)1/2 ^ (32) 

771 —— 2 

where in the last step we have used Eq. (l30l) . 

From the above we find 

02(0^)'50/«|O) = 2Q{6R + 2Q + 3){elf-\el)^\0) 

+ 3i7(i7-l) (^t)Q+2(Qt)P-2|o)_ (33) 
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Note that this implies, for the special case R = 0, 

0502(05)‘5|O) = 2Q{2Q + 3)(0j)<3|O) , 


The state dJTl i can be expressed in terms of the basis 
|n2, ni, no, n_i, n_2), where n^ is the number of particles 
(34) in the state m. We get 


a result which we shall see again in Sec|IV]of this SM. 

We can now derive the exact manybody wavefunction for 
the region C when TV is a multiple of 3 . The state |4'3fl) with 
N = 3R particles and t = N with N = 3 i? being a mul¬ 
tiple of 3 {i.e. 0202 |4'3n) = 0 hence A = 0 ) can then be 
constructed as 

\^zr) = b^{Ql)^ + + ... 

+ 6fc(0j)3'=(0t)«-2fc + ..., (35) 


where we have 


|4'2q) = 


3 X 2<3Q! 


1 1/2 


(2Q + 3)!! 

X |fc2,fcl,2fco,fci,fc2), 


Q 

ko,ki,k2 


[(2fco)!]^/^ 

2^-0 fco! 


(42) 


where the sum is over all non-negative integers kg, ki, k 2 with 
the restriction (denoted by the prime) fco + + ^2 = Q- 

The density matrices are obtained by operating on 

|T' 2 q) and then evaluating the appropriate inner products. The 
required sums are evaluated below. 

We show here how to evaluate the sums involved. They are 
of the form 


bk+i ^ {R-2k){R-2k-l) 

bk l2{k + l)[2R-2k-iy ^ ’ 

We also note here that since [ 02 , N{N + 3) — 0202 ] = 0, 
the states (02)^|4'3i{) have the same quantum number r = 
3i? though different particle numbers N = 2Q -b 3i?. From 
these we obtain the exact many-body ground states in region 
C of Fig. 1 in the main text. 


IV. REDUCED DENSITY MATRIX CALCULATION 
We here consider the N = 2Q singlet state 

|v1/2q) = ^^0tQ|O), (37) 

where g{Q) is a normalization constant. g{Q) can be evalu¬ 
ated (see also below) by the repeated use of the commutation 
relation [0, 0^] = 2(2TV-|-5) where N is the number operator. 
We then obtain 


5 (g) = 2^Q!(2Q + 3)!!/3. (38) 

Some special values are; 5 ( 1 ) = 10, g{2) = 2^ x 5 x 7, and 

0 ( 3 ) = 2'^ X 3^ X 5 X 7. It turns out that g{Q) = /(I, Q) of 

The expectation values for the two-particle density ma¬ 
trices {^ 2 q\^jm^Jm\^ 2 q) needed can be read off from 
the energy i? = ^ [aN{N — 1) -|- j3{F{F -b 1) — QN) -b 7A] 
where A = 2Q{2Q -b 3) since this must also be i? = 
iZ]jm5f(4'2q|4jm"^Jm|4'2q) where the sum over J is 
for 0, 2,4 only. We have 

(T' 2 Q|iSo^o|^ 2 Q) = 2g(2g + 3)/5 (39) 

{^2q\AI^A 2 m|'I' 2 q) = 8 Q(Q — l)/35 (40) 

(^ 2 Q|ilM A4 m|4'2q) = 8 g(g — l)/35 (41) 

The equality between the values between J = 2 and J = 4 is 
due to the special properties of the state | 4 ' 2 q). In below we 
also show an alternate derivation of Eq. (I39ll411 i 


Q 

SQ = J2h, 

k=0 


(43) 


where bk are the products of polynomials in k with Ck = 
22 k^k\y ~ ^^ 2 ^fcV'' • notice that if the function 

fiy) = known, then (by straight-forward ver¬ 

ification) Sq is simply the coefficient of of the function 
F{y) = /( 5 )/(l-y). Now, we note that /i(y) = ^^y^ 

is given simply by (1 — Hence the sum Si^q = 

Y^k=o given by the coefficient of (1 — y)~^^^, and 
hence 


Q 

k—0 


{2Q + iy.\ 
2QQ\ 


(44) 


Similarly, for the sum S 2 ,q = Y^k=o function 

h{y) = can be obtained from y-^fi{y) = 

|(1 — Hence S 2 .Q is the y^ coefficient of ^ 2 ( 2 /) = 

|(1 — y)~^^^, and hence 


Q 

S2,q = ^ kck 

k=0 


Q{2Q + iy.\ 
3 2'7g! 


We can proceed similarly to get 


g(g- 1 ) ( 2 g + i)!! 

5 2QQ\ ’ 


Q 

S 3 ,Q = ^k{k- l)ck = 

k=0 

Q 

Si,Q = Y,k[k-l){k-2)ck, 

k=0 

Q(Q-l)(g- 2 ) ( 2 Q + 1 )!! 

7 2'2g! 

With the above sums, we can also obtain 


y~^(g +1 - k)ck 




^(g - fc)(g -k + i)ck 


k=0 


1 (2g + 3)!! 

3 2<3g! ’ 

4g (2g + 3)!! 
15 2Qg! ’ 


(45) 


(46) 


(47) 


(48) 

(49) 
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(using {Q - k){Q - k + 1) = k{k - 1) - 2Qk + Q{Q + 1)), 
and 


Q 

- k - 1)(Q - k){Q - k + l)cfe 

k=0 

8Q{Q-1) (2g + 3)!! 

35 2<3Q! ’ 


whereas 


(^2q|44“i®iI^2q) = 0. (55) 

The last result can be most easily seen when T' 2 Q is expanded 
in the number basis. This is reflected in the equality between 
Eqs. drol l and (HTl i. 

For M=l, there are two operators aiuo and d 2 d_i. We have 


(using {Q — k — l)(g — k){Q — k + 1) = —k{k — l){k — 
2) + 3(g - l)k{k - 1) - 3Q{Q - l)k + {Q- 1)Q{Q + 1)). 

We demonstrate the use of the above relations by check¬ 
ing here the normalization of |' 1 ' 2 q)- (^' 2 q|'I' 2 q) is given by 


3x2'3Q! 

(2Q-I-3)!! 


y^q 

Z^fco 


k^ ,k- 


'cfc„. Due to the restriction ko + ki + 


k2 = Q, the sum is therefore given by I]fco=oEfei=o" 

= Ck{Q - k + 1), which is from Eq. (gOll. The 

density matrices are obtained by first operating a™ or dmi am^ 
on | 4 ' 2 q) and then evaluating the appropriate inter-products 
with the help of the above formulas. 

We list here also the two-particle density matrices 
{'^ 2 Q\aln^aln^am 3 arrn\^ 2 Q)- We list them starting from the 
largest M = mi + m 2 = m 3 -I- mi. For M = 4, 


(T'2Q|a|4a2d2|4'2Q) = 8Q{Q - l)/35. (51) 


For M = 3 


(T'2Q|a^aIaia2|4'2Q) = 4Q(Q - l)/35, (52) 


corresponding to (^' 2 q|^ 1 m^ 4 m|^' 2 q) = 8Q{Q - l)/35. 
ForM=2, we have two operators 000,2 and oidi and their con¬ 
jugates. We obtain 

(4-20144^002|«'2q) = 4g(Q - l)/35, (53) 

(4'2Q|aIaIaiai|T'2Q) = 8 g(Q- l)/35, (54) 

I 


(T'2Q|a|aJdodi|«'2Q) = {'H 2Q\a\a^-ia-ia2\'i! 2 q) : (56) 

= 4g(Q-l)/35, (57) 

and there are no cross elements. Similar remarks we made for 
the M = 2 sector also applies here. 

For M=0, there are three operators, 020,-2 and oio-i 
and dodo. This part of the density matrix is, with rows and 
columns in order of these three operators, given by 

on ( (4Q + 3) -(2Q + 5) (2Q + 5) \ 

^ -(2Q + 5) (4Q + 3) -(2Q + 5) . (58) 

\ (2Q + 5) -(2Q + 5) (6Q + 1) / 

The equality between the first two diagonal elements, as well 
as among the off-diagonal elements except signs, follows 
from the fact that |4'2 q) is invariant under d ±2 a±i up 
to a sign. From the above formulas, we can recover Eqs. ( l39l l- 

(EB. 

It is straightforward to obtain the number fluctuations from 
above. {N 2 N 2 ) = {NiNi) = {N 2 N_ 2 ) = {NiN_i) = 
2Q{AQ + 3)/35, (iVoiVo) = 4Q(3Q + 4)/35, {N 2 N 1 ) = 
{N 2 N 0 ) = {N 2 N- 1 ) = (NiNo) = 4Q(g - l)/35. The 
above expressions are valid also when replacing m by — m. 
{N 2 N 2 ) = {N 2 N- 2 } etc follows immediately also from (l42l) . 
since for each of the states on the right-hand-side, 712 = ^ 1 - 2 - 
We now consider the density matrices for |4'p)aD in large 
N limit. We have 


(UmiUmiUma 07714 ) ~ N(N 1) 




I N-2 




1 N 


(59) 


(We leave out the explicit labels |4'p)a„ to simplify our no¬ 
tations). For large N, the overlap [J2m ‘P*mi^l)‘Pm{^2)]^ is 
negligible unless fli is very close to 1 ) 2 - (rigorously speaking, 
also r 22 — but we can check easily that this does not affect 
the following argument). Hence we can identify the H’s in the 
arguments of ipmi in the integrand of the numerator. 

Canceling the common factors (the normalization coefficient 


is Ok N ^ for large N) we are left with 

(Umi Umi Urns ^7774 ) 

= 7V(7v-i) (0)^27,73(60) 

Jn 

Now we observe that, for our state, ip^{Ct) = (—l)'"<p_m(^), 
since = (-l)”"+"d )^^77 where /3, 7 ) = 

g- 7 ( 777 'a-i- 7777 )^^) Therefore the integral above is the 
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same as 

Jn 

We notice that this latter integral is the same as the one that oc¬ 
curs in our evaluation of the coefficient of ®-m, 2 ®rn, 3 ®m 4 
for the wavefunction |'l'p)a« for four particles, except combi¬ 
natorial factors. For example, the value of (a 2 a|aia 2 ) is just 
—N{N — 1) times the coefficient of a^_ 2 a^_ia\a\ in |4'p)a« 


(see Eq. (flSl l) divided by 4!. Using our previous calcula¬ 
tions we therefore obtain, in the large N limit, (d 2 d|dia 2 ) = 
iV^/35, (a 2 aJaoa 2 ) = iV^/35 etc. It is again most economi¬ 
cal to express the final results using Ajm- We get (AJqAqo) = 
iV^/5 and {A\j^A 2 m) = {A\j^AiM) = 2iV^/35. Hence the 
terms in two-particle density matrix in the state obtained 
by angular average is the same as that of Eq. dJTl i. and the 
differences arise only in lower powers in N. Therefore the in¬ 
teraction energies per particle for these states are equal except 
for terms that are of order 1. 


